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Numerical study of staggered fermion on anisotropic lattices* 

Kouji Nomura^, Takashi Umeda'', and Hideo Matsufuru*^ , 
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We study calibration procedures of the staggered quarks on anisotropic lattices in the quenched approximation 
and in Nf = 2 dynamical simulations. For the calibration conditions we adopt the hadronic radii ro and the 
meson masses in the temporal and spatial directions. On the quenched lattice, we calibrate the quark field and 
compare the result with the result determined using the meson dispersion relation. In dynamical simulations, we 
determine the anisotropy parameters 7g and 7f simultaneously within 1% accuracy at renormalized anisotropy 
^ = Os/at = 4. 



1. Introduction 

Anisotropic lattices realize small temporal lat- 
tice spacings while keeping modest computational 
costs. The technique is useful in various fields 
of lattice QCD simulation: At finite temperature 
the large number of degrees of freedom in Eu- 
clidean time direction leads large number of Mat- 
subara frequencies, which is efficient for calcula- 
tion of the equation of state ^ and for analy- 
ses of temporal correlation functions |2I3I4| . The 
large temporal cutoff is important for relativistic 
formulations of heavy quarks on lattices |5] . It is 
also convenient for the correlators in which noises 
quickly grow as time such as glueballs |4I6| and 
negative parity baryons 0- 

However, uncertainties of anisotropy param- 
eters bring additional errors into the observed 
quantities. For precise calculations, we need to 
tune anisotropy parameters with good statistical 
accuracies and to control the systematic errors 
in the continuum extrapolations. In this work 
we discuss the tuning procedures of anisotropy 
parameters for the staggered fermions in the 
quenched and dynamical QCD simulations. Our 
goal is to tune anisotropy parameters within 1% 
statistical errors at each set of (3 and mo- In the 
following, we focus on lattices with the renormal- 
ized anisotropy = as/at — 4,. 



2. Lattice action 

In the numerical simulations, we adopt Wilson 
gauge and staggered quark actions on anisotropic 
lattices. The Staggered quarks have several ad- 
vantages over the Wilson-type quarks in studies 
related to the chiral symmetry and in simulations 
at lighter quark masses. The gauge action is 
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where 7g is the bare gauge anisotropy parameter. 
The quark action is defined as 

Sf =Y.^{x)K{x,y)x{y). (2) 
K{x,y) = 6x,y 

-~i^<T Vt.{x) [Ur{x)S^ y_^^ -Ujix- i)6^ y^i'^ (3) 

i 

where is a bare fermionic anisotropy parame- 
ter, rifj,{x) the staggered phase, = l/2mo with 
mo the bare quark mass in spatial lattice units. 

Although these anisotropic lattice actions have 
been adopted in Ref.^, they did not discuss the 
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systematic uncertainties due to the anisotropy to 
the sufficient level for precision computations. 

3. Calibration procedures 

The anisotropy parameters 7g and 7i? should 
be tuned so that the physical isotropy condition, 
Cg(7g'7f) = ^Fhh'^*F) = ^ holds at each set 
of (3 and mo, where ^gItg, 7-f) and Cf(7g,7f) 
are the renormalized anisotropics defined through 
the gauge and quark observables, respectively. In 
quenched simulations, one can firstly determine 
the 7G independently of jf, and then tune 
for fixed jq. On the other hand, in dynamical 
simulations, 7g and ^f need to be tuned simul- 
taneously. 

3.1. Gauge sector 

For the calibration in the gauge sector, we de- 
fine Cg(7G)7f) through the hadronic radii rp 
in the spatial and temporal directions. Since we 
set the lattice scale via ro , the physical anisotropy 
is kept in the continuum extrapolation. This is 
an advantage to avoid that the systematic errors 
brought by the anisotropic lattice remain in the 
continuum limit. 

The value of rg is defined with the static q- 
q potential through the relation r§i^(ro) — 1.65, 
where F{r) is the force between quarks. Then 
is defined with the ratio of tq's in the temporal 
and spatial directions: ^g 
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In quenched simulations, we can apply the 
Liischer-Weisz noise reduction technique |1U| and 
determine £_g at 0.2% level of accuracy [T^ . 
In dynamical simulations, we instead apply the 
smearing technique to the anisotropic three- 
dimensional planes which enables determination 
of ^G at 1% level of accurately. 

3.2. Quark sector 

In the fermion sector, we define £,f with the 
ratio of meson masses in temporal and spatial di- 
rections: ^F = ITT- /'ITT' 12]' We use the pseu- 
doscalar channel, since in that channel the os- 
cillating components of the correlators disappear 
quickly and hence statistical errors of masses are 
small. 

Alternative definition of the £,f makes use of 
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Figure 1. The result of calibration of the stag- 
gered quark field in the quenched simulation. 



the meson dispersion relation 



(4) 



Here E(jp) and m are in temporal lattice units and 
p is in spatial lattice units, and hence appears. 
Pi = l-KTLilLi (i =1,2,3), where Li is lattice size 
in i-th direction. We use n = (0, 0, 0), (0, 0, 1) to 
define ^i?. 

At present stage of this work, we adopt the 
former procedure, with the ratio of the temporal 
and spatial meson masses, as the main calibration 
procedure for the quark field, because of smaller 
statistical uncertainties than the latter definition. 
In quenched simulation, we compare the results 
with these two procedures. 

4. Numerical results of calibration 

4.1. Quenched simulation 

The simulation is performed on a 12^ x 24 x 
96 lattice a.t /3 — 5.75 and 7g — 3.136 which 
correspond to the spatial cutoff a^^ —1.1 GeV 
and the renormalized anisotropy ^g = 4 T2, . The 
statistics are 224 configurations. The bare quark 
mass range is from 0.01 to 0.5 which corresponds 
to the range mps/mv from 0.33 to 0.90. 

We calculate the meson mass ratio at several 
values of 7^ and interpolate them to 7^^ at which 
^F = 4 holds. The result is displayed in Fig. ^ 
For our lightest quark mass case, is deter- 
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Figure 2. The result of calibration in the dy- 
namical simulation at P = 5.3, mo — 0.1. The 
top and bottom panel show the dependences of 
7g on and on ^p, respectively. 

mined within 0.5% statistical error. 

In Fig. m we also show the result of calibra- 
tion using the meson dispersion relation. The re- 
sults of two procedures are consistent only in vary 
light quark mass region. To understand the dis- 
crepancies of two procedures and large change of 
7^ from the calibration with mass ratio at larger 
quark mass region are subjects of future work. 

4.2. dynamical simulation 

The simulation is performed with R-algorithm 
with Nf = 2 on lattices of the size 10^ x 20 x 80 
at 4 values of (3 from 5.3 to 5.45 and the quark 
mass mo = 0.1. The statistics are from 300 to 800 
trajectories after 200 for thermalization. In the 



following we present only the result at /? = 5.3 on 
which the spatial cutoff is = 0.69 GeV and 
mps/mv — 0.6. 

We determine 7^ and 7^, by a linear interpo- 
lation to ^ = 4 using the 6 parameter sets of 
(7G:7f) assuming the following relations: 
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ao,ai,a2j = 

[&0,6l,52] = 

We find that 



The result of calibration with the meson mass ra- 
tio is presented in FigI21 

The fit results in 

[4.000(24), 0.702(92), 0.075(92)], 
[4.000(28), 0.220(08), 1.225(08)]. 
(^F has small dependence on 70 and on ^f- 
By solving Eq. 10), we obtain 7^ = 2.885(35), 
1p = 2.955(31). We can tune the parameters 70 
and 7f with the 1% level of statistical errors. 

This calculation was done on a Hitachi SR8000 
at KEK. H.M. and T.U. thank JSPS for Young 
Scientists for financial support. 
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